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Ab initio calculations have been performed on hexagonal layers of M2O3 (M being several tran¬ 
sition metals of the 5 d series) sandwiched by a band insulator such as AI2O3 that provides the 
honeycomb lattice where the 5 d electrons reside. This corundum-structure-based superlattice is the 
most obvious way to design a honeycomb lattice with transition metal cations avoiding the use of 
largely polar surfaces. We obtain that this system supports the presence of Dirac cones at the Fermi 
level that open up with the introduction of spin-orbit coupling at various fillings of the 5 d band. The 
DFT calculations performed in this work show that the 5 d 5 situation is always a trivial insulator, 
whereas the 5 d s filling presents topological insulating configurations which evolve into a trivial state 
with increasing tensile strain or on-site Coulomb potential U. However, LDA+U calculations show a 
stable antiferromagnetic solution for the 5d 8 case at every U value, which would break time reversal 
symmetry and could affect the topological properties of the system. We also discuss the similarities 
with the buckled honeycomb lattice obtained using perovskite (111) bilayers previously studied in 
literature, in particular for the 5d 5 and 5 d s configurations. This work provides some clues on the 
stability of topological phases using metal oxides in general. 


I. INTRODUCTION 


In recent times, the interest in topological 
insulators^ has become ubiquitous in condensed 
matter physics literatureThe first observed 
examples were found in the classic thermoelec¬ 
tric materials^ and their nanostructures^- 7 but 
soon the list expanded to include oxides A— com¬ 
pounds based on Bi, Te and Sb^ - — and even 
Kondo-type insulators —r— Various types of unique 
topological phases have been explored, includ¬ 
ing the interface between a topological insulator 
and a superconductor topological ordering in 
superconductors^ the existence of so-called topo¬ 
logical Mott insulators topological crystalline 
insulatorsand the list continues. We are on 
the way to start thinking about building devices 
such as field effect transistors out of topological 
insulators ; 35 ' 36 and yet their properties are still in¬ 
triguing and good candidates for applications in 
electronics ( e.g. having a large band gap2£) or in 
quantum computing (e.g. by producing Majorana 
fermion s 38 ' 39 ) still are to be found experimentally. 
Topological insulators (TI) constitute a new state of 
matter that has forced us to rewrite our basic solid 
state physics textbooks. They are insulating mate¬ 
rials in their bulk state but at the same time they 


present dissipationless chiral channels on the surface. 
Those edge states present an absence of back scatter¬ 
ing that opens new possibilities for charge and spin 
transport in nanostructures. They are also topolog¬ 
ically protected by certain symmetries of the one- 
particle Hamiltonian that describes the properties 
of the system. The quest in under way to find new 
topological insulators that could be stable and us¬ 
able. One possible path is to try to design them 
based on oxides. These have the advantage of being 
very well known materials, stable in ambient con¬ 
ditions and can be grown in various nanostructures, 
such as thin films or multilayers using standard tech¬ 
niques. One such route has been explored recently 
by various authorsi 10 ' 12 ' 40- — In particular, creating 
a honeycomb lattice using a perovskite bilayer— has 
produced interesting predictions for various 5 d tran¬ 
sition metal compounds being the active part of the 
multilayered system. We propose here an analo¬ 
gous design but building the honeycomb lattice in 
a more natural way, avoiding the use of a largely 
unstable surface of the perovskite, the highly po¬ 
lar (111) surface that only recently has started to 
become accessible for experimentalists— We will 
use instead a non-polar surface based on the corun¬ 
dum structured that provides intrisically a transi¬ 
tion metal-based honeycomb lattice. 

The paper is organized as follows: in Section |TT] 




FIG. 1. (Color online) (a) and (c) Scheme of the corundum structure multilayer (M 203 )/(Al 203 )s studied where the 
active single layer M2O3 is highlighted, (b) Top view of the honeycomb lattice formed by the 5 d TM cations M. (d) 
Reciprocal lattice showing the path of the band structures presented. 


we will present the details of the ab initio calcu¬ 
lations, in Section eh we will describe the system 
under study, in Section IIVI we will analyze the re¬ 
sults of our calculations for the most interesting 
compounds, suggesting the analogies and differences 
with results already available in literature: Section 
IIV Al describes the 5 d 5 and 5 d 8 systems as possible 
candidate band structures to host non-trivial topo¬ 
logical properties, in Section IIVBI we provide a de¬ 
scription of the magnetic properties of these systems, 
in Section HV Cl we provide a detailed study of each 
case as a function of strain and in Section llV Dl we in¬ 
clude the effects of correlations at the LDA+U level. 
Finally, in Section[V]we will summarise the main im¬ 
plications of our work. 


II. COMPUTATIONAL DETAILS 

Ab initio electronic structure calculations based 
on the density functional theor y 49 ! 50 (DFT) have 
been performed using the all-electron full potential 
WIEN2I-—> software. The calculations were per¬ 
formed for different off-plane and in-place lattice 
parameters of the corundum structure, all of them 
were relaxed using the full symmetry of the original 
cell. The exchange-correlation term depends on the 


case. Uncorrelated calculations were performed us¬ 
ing the Generalized Gradient Approximation (GGA) 
with the Perdew-Burke-Ernzerhof schemed The Lo¬ 
cal Density Approximation+U (LDA+U) in the so- 
called “fully localized limit”— was utilised in order 
to explore the effects of on-site Coulomb repulsion in 
the 5 d manifold. For selected configurations, GGA£ 2 
and LDA-^ calculations are shown to yield qualita¬ 
tively equivalent results. This will be presented and 
discussed in Section Hvl 

The calculations were performed with a k-mesh 
of 8 x 8 x 8 and RmtK-max = 7-0 (the product of 
the smallest muffin-tin radius and the maximum k- 
vector of the plane-wave expansion determines the 
size of our basis set). Spin-orbit coupling (SOC) 
was introduced in a second variational manner by 
using the scalar relativistic approximation^ these 
calculations were performed with a modified k-mesh 
of 8 x 8 x 2. In the GGA calculations of the bulk 
structures of M2O3 of Os and Au, the muffin-tin 
radii R mt (in a.u.) were set to 1.92 for the metal M 
and 1.57 for the O. In the (Au 2 03)/(A1 2 03)5 mul¬ 
tilayer, the R mt values used were: 1.78 for Au, 1.68 
for Al and 1.5 for O; for the (0s 2 03)/(Al 2 03)s mul¬ 
tilayer: 1.78 for Os, 1.68 for Al and 1.58 for O; in the 
(Ir 2 0 3 )/(A1 2 0 3 )5 multilayer: 1.92 for Ir, 1.61 for Al 
and 1.61 for O. 
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III. THE SYSTEM 


0.2 


We can see in Fig. (T] the multilayered struc¬ 
ture of the system under study, where a single layer 
of M 2 O 3 , M being a 5 d transition metal (TM), is 
sandwiched by a wide gap band insulator, such as 
with the so-called corundum structure 
that provides the symmetry of the multilayer. In 
that geometry (Fig. |T] (b)), the TM cations sit 
on a honeycomb lattice, resembling the situation 
in graphene or in other oxide-based heterostruc¬ 
tures, such as the aforementioned perovskite bilayers 
grown along the ( 111 ) direction, but without the ad¬ 
ditional bucklings caused by the bilayer structure, it 
is a single-layer here that provides the honeycomb. 
For that reason, this system is a more natural real¬ 
ization of a 5 d metal-based honeycomb lattice util¬ 
ising non-polar oxides for its growth. If a similar 
electronic structure occurs, then comparisons can be 
drawn with other similar systems that sustain topo¬ 
logically non-trivial phases. We provide a descrip¬ 
tion of these in the following Section. 


IV. RESULTS 


A. General overview of the electronic structure 


DFT calculations were performed for various 5 d 
TM cations on the M 2 O 3 sandwiched layer. We can 
see in Fig. [2] the band structure of the multilayer 
around the Fermi level with and without spin-orbit 
coupling for various fillings of interest: half-filled e g 
levels (AU 2 O 3 , panels (a) and (b)), one hole in the 
t 2 g multiplet (OS 2 O 3 , panels (c) and (d)) and 5c ? 6 
configuration (I 1 - 2 O 3 , panels (e) and (f)) that has a 
full t 2 g shell but shows a Dirac cone at 5 d, 5 filling. 
We have studied all the other possible fillings but 
none of them show anything close to a Dirac point 
near the Fermi level. 

In the perovskite-based bilayer system, two of 
these fillings (5 d 8 and 5 d 5 ) have been predicted to be 
topologically non-triviaL 10 ’ 12 Here we see a similar 
band structure. However, for the 5c ? 8 case, a band 
crossing at the Fermi level at the M point is removed 
by SOC, whereas in the perovskite bilayer case it oc¬ 
curs at K j 10 i 12 This suggests that the tight binding 
model for this unit cell is different. This is consis¬ 
tent with the different hopping paths in the corun- 
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FIG. 2. (a) and (b) show the band structure ob¬ 

tained within GGA for the case without (left)/with(right 
panel) spin-orbit coupling, respectively, for the 5 d s filling 
(AU 2 O 3 ) at a = dAi 2 o 3 and optimized c lattice param¬ 
eter c = 26.51 A. Only the four e g bands are shown, 
(c) and (d) show similar calculations for the 5d 5 filling 
(OS 2 O 3 ) with c = 25.78 A, where only the t 2 9 states are 
shown, (e) and (f) show the t^g manifold for the 5 d 6 
filling (R 2 O 3 ) with c = 25.50 A, where the Dirac cone 
opens up with SOC at the 5 d 5 filling. 


dum unit cell in comparison with the perovskite case. 
Yet, it is worth analysing the similarities in terms of 
their topological properties. 

Also, for the 5c ? 5 case, we notice that the bands are 
considerably broader and hence the relative strength 
of SOC will be smaller compared to the perovskite- 
bilayer case. We will explore below what hap¬ 
pens when the SOC strength is artificially enhanced. 
Moreover, for the Os-based multilayer we see the 
gap is not affected by SOC. However, it looks like 
the system can be close to a topological transition 
via increasing SOC, reducing hopping parameter f, 
applying strain, etc. since the Ir-based multilayer 
shows a Dirac cone at 5c ? 5 filling which responds to 
SOC by opening a gap as shown in Fig. [ 2 ] (e) and 
(f). We will explore further all these possibilities 
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FIG. 3. (Color online) (a) Representation of the magnetic moment evolution with U, both for 5 d 5 (black circles) and 
5d 8 (red squares) fillings. The shaded region highlights the magnetic moment calculations inside the typical realistic 
values of U for 5 d electrons. Notice the similarity in the evolution of the magnetic moments between the two systems 
studied, but just at two different U regimes, (b) Schematic representation of the AF solution in the honeycomb plane. 


throughout the paper. 

Due to the presence of time reversal symmetry 
(TRS) and inversion symmetry (IS), the topologi¬ 
cal invariants of the band structure can be easily 
computed in the bulk state and we can determine 
whether the system is a TI or a trivial insulator by 
following the simple prescription of Fu and Kane. 
This computes the Z 2 invariant of the system utiliz¬ 
ing only the single-particle eigenstates at the time- 
reversal-symmetric k-points of the Brillouin zone^ 
as long as TRS and IS are retained, the latter oc¬ 
curring for all the calculations presented. 

By means of this rule, we will be able to analyse in 
detail whether this multilayered system (when it is 
TR symmetric) is topologically trivial or non-trivial, 
depending on the Z 2 invariant that results from the 
calculations. 


B. Importance of magnetism and correlations 

Before we apply the so-called Fu-Kane^, which is 
only valid when TRS and IS are present, let us recall 
that this requieres that the system must have a non¬ 
magnetic solution, since IS is always retained in our 
calculations. We should clarify in what regime these 
multilayers will present a non-magnetic solution as 
most stable. For that sake, we have studied the evo¬ 
lution of magnetism with the on-site Coulomb re¬ 
pulsion U by means of the LDA+U method for the 


(M 2 O 3 ) /(Al 2 03)5 system by setting up a nearest 
neighbor antiferromagnetic (AF) coupling between 
the transition metals in the honeycomb plane (Fig. 
OH (b)). Being rather delocalized, the 5d electrons 
in these systems are usually treated with values of 
the on-site Coulomb repulsion on the order of 1 — 3 
eV, the correct values being difficult to estimate, we 
present results calculated at various U values. 

The evolution of magnetism with the Coulomb in¬ 
teraction is summarised in Fig. [3] for the two fillings 
we will focus along the paper in. We see that as 
U increases, an AF solution becomes stable for both 
compounds, the onset point being at different values 
of U. 

The Os-based system presents an AF ordering for 
U > 3.5 eV (a sizable value for these 5d-electron sys¬ 
tems) that would only enhance the gap without pro¬ 
ducing any crossover to a phase with different topo¬ 
logical properties. In contrast, the Au-based system 
has a stable AF solution for any U > 0. Even though 
our DFT calculations predict an AF phase to be sta¬ 
ble, let us point out that AU 2 O 3 in the bulk form 
is itself non-magnetic42 The same non-magnetic re¬ 
sult is obtained from our own calculations for bulk 
AU 2 O 3 with corundum structure. Also, another 
Au 3+ compound like LaAuC >3 was found experimen¬ 
tally to be diamagneticThus, it is interesting here 
to study the topological properties of a non-magnetic 
phase as well, as we will carry out in detail below. 
Yet, more studies are required to understand if a 
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hexagonal single-layer AU 2 O 3 would sustain mag¬ 
netism or not. In case it does, such system would 
be very interesting to study the effects of strong cor¬ 
relations in those very narrow bands produced by 
spatial confinement. Previous works based on per- 
ovskites bilayer s 12 ! 59 have also studied the magnetic 
solutions of the 5 d 5 and 5 d 8 compounds. For the 
5 d 5 filling case ^ 2 an AF solution becomes stable 
for U > 3 eV breaking TR and IS which is simi¬ 
lar to the results obtained for the corundum-based 
system presented here. On the other hand, the 5 d 8 
system^ presents a narrow window inside the AF 
region (U > 2.1 eV) where the system remains in a 
topological insulating state. By increasing the pa¬ 
rameter U, this phase is destroyed and the system 
evolves (through gap closing) into a trivial insulator 
with AF ordering. In this corundum-based system, 
the e g bandwidths are substantially reduced with 
respect to the perovskite-based case, promoting the 
strength of the magnetic coupling. Yet, in principle 
magnetism need not break the topological insulating 
phase. Ref. [l 2 j shows that breaking inversion sym¬ 
metry does not immediately destroy the topological 
insulating phase in the case of perovskite bilayers 
either. The actual topological features of the AF 
phase for these corundum-based systems are beyond 
the scope of this paper. 


C. Detailed analysis case by case. 

Once we have clarified what DFT calculations pre¬ 
dict about the stability of magnetism in these sys¬ 
tems, we focus on analyzing the topological proper¬ 
ties of non-magnetic phases and their strain depen¬ 
dence carrying out calculations at different values of 
U but retaining a non-magnetic solution where TRS 
is not broken. 


1. (Au2 03)/(Al 2 03)s system 

Let us first discuss the properties for the 5 d 8 fill¬ 
ing. This is a half-filled e g shell, with all the t 2 g 
levels completely occupied and far below the Fermi 
level. We will forget about them in our description. 

Figured] (a) shows the evolution of the gap and 
the total energy of the system (Fig. [4] (b)) with 
in-plane strain (lattice parameter a). The picture 


shows how, as a increases, the topological gap de¬ 
creases to zero and reopens again as a trivial insula¬ 
tor. The total energy curve shows that the crossing 
point (gap zero) would be very close to the equi¬ 
librium structure. Due to the lattice size mismatch 
between AI 2 O 3 and AU 2 O 3 (which is quite substan¬ 
tial; a GGA-based structural optimization in bulk 
AU 2 O 3 with corundum structure yields a = 5.48 A 
and c = 15.6 A, that is a sizable in-plane mismatch 
of about 15 %), the lattice prefers to deform its vol¬ 
ume increasing the a lattice parameter compared to 
that of AI 2 O 3 . If these multilayers can be grown epi¬ 
taxially retaining the in-plane lattice parameter of 
the substrate (AI 2 O 3 , a = 4.75 A), the system would 
be well inside the topologically non-trivial region of 
the phase diagram. However, our total energy calcu¬ 
lations suggest a more stable situation would occur 
closer to the transition point with a much smaller 
band gap, and also that epitaxial growth on a sub¬ 
strate with larger a could tend to reduce the topo¬ 
logical gap. In any case, tensile strain drives the 
system towards a trivial insulating phase whereas 
compressive strain induces a topological insulating 
phase. 

Figured] (c)-(e) shows the evolution of the band 
structure with strain. We see there the band struc¬ 
ture for three different values of the in-plane a lattice 
parameter, panel (c) with the a of AI 2 O 3 and pan¬ 
els (d) and (e) for two larger a values. For each 
a parameter, the c-lattice parameter was optimised. 
We see that the gap at the M point closes as a in¬ 
creases and then reopens again. When it does so, 
the v value of the bands around the Fermi level has 
changed. By computing the Z 2 invariant following 
the Fu-Kane criterion, we see that the system is a 
topological insulator for a = 4.75 A, the gap be¬ 
comes zero at around a = 4.89 A, and it reopens as 
a trivial insulator at larger a values. 


2 . ( Os 2 O3)/[Al 2 03)5 system 

A similar analysis has been carried out for the 
(OS 2 O 3 )/(Al 2 03)5 case. In this 5 d 5 filling (one hole 
in the t 2 g manifold), corresponding to the nominal 
occupation of the Os 3+ cation, the existence of a 
topological insulating phase is not observed. All the 
configurations studied were found to be trivial in¬ 
sulators (or metallic phases). We have optimized 
the atomic positions and the volume (a and c lat- 
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FIG. 4. (Color online) Evolution of the band gap (a) and the total energy (b) for the (Au 203 )/(Al 203 )s system. 
All the GGA calculations were performed at the optimized lattice parameter c with relaxed atomic positions. The 
inset zooms at the region close to the energy minimum. Band structure evolution (only the e g manifold is shown) 
for the 5d 8 filling with strain: (c) aAi 2 o 3 = 4.76 A, (d) clgs = 4.89 A and (e) a = 4.95 A as the systems evolves 
from a topological to a trivial insulator by closing and reopening the gap. The topological invariants of each band 
are calculated and shown. 


tice parameters) and the ground state was found for 
a = 4.80 A, very close to the value for AI 2 O 3 . A wide 
range of both tensile and compressive strains were 
tested by varying the in-plane a lattice parameter 
(optimizing c in every case). None of the strains con¬ 
sidered yielded a topological insulating phase. How¬ 
ever, we observed that a Dirac cone at the K point 
near the Fermi level that opens up with SOC ex¬ 
ists for the 5 d 5 filling when the active layer is H 2 O 3 , 
nominally a 5 d 6 state (if the chemical potential were 
shifted to one electron less). This could indicate that 
the system is close to a transition towards a topo¬ 
logically non-trivial solution, which we have tried 
to explore. This 5 d 5 filling has been predicted to 
host topologically non-trivial phases for the buckled 
honeycomb lattice constructed using a perovskite bi¬ 
layer along the (111) direction. We can see that the 
situation changes here for OS 2 O 3 as an active layer, 
whose band structure including spin-orbit coupling 
can be seen in Fig. [5] (a). This shows only the 
t 2 g states, with one unoccupied band just above the 


Fermi level. The band width of the t 2 g multiplet is 
quite large, of about 3 eV, substantially larger than 
in the perovskite bilayer system that hosts topolog¬ 
ically non-trivial phases (where the band width is 
about 1.5 eV)i 10 i 12 In the corundum case, the hop¬ 
pings occur all of them in one plane (compared to 
the buckled honeycomb lattice in the case of the per¬ 
ovskite bilayers) and broader bands show up. In 
both cases, the top of the valence band is at the K 
point. The fact that the ratio A so /t (being A so the 
SOC strength parameter) is small implies that SOC 
is not strong enough to promote a band inversion 
that could lead to a topologically non trivial phase. 
Still, two-dimensional confinement leads to an insu¬ 
lating phase which survives for a significant range of 
strain effects. 

In order to test how far this system is from a 
spin-orbit induced phase transition we will enhance 
the strength of the SOC and calculate the electronic 
structure for various values of A so , where the stan¬ 
dard spin-orbit term in the Hamiltonian was consid- 
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FIG. 5. (Color online) (a) Band structure calculated with spin-orbit coupling showing the i 2 9 manifold of the system 
(OS 2 O 3 )/(AI 2 O 3)5 for the ground state configuration with a = 4.80 A and c = 13.64 A. The topological invariants 
calculated show a trivial insulating state that holds for all the lattice parameters studied. Observe that the gap is 
sizable (about 200 meV) in this case, too big to be easily reduced by applying strain, (b) Evolution of gap with the 
SOC strength parameter A so /Ao (the artificially modified value divided by the actual one) at k-points M and K. The 
gap becomes zero at approximately A so /Ao « 2.1 where the gap presents a discontinuity at both k-points studied. 
Several calculations with slightly different values of A so /Ao around 2.1 were performed but no intermediate value of 
the gap in the interval was observed, there is a discontinuous change there at both k-points simultaneously. Even 
though there is a drastic phase change observed, all the calculations preformed show a trivial insulating phase for 
(0S20 3 )/(A1 2 0 3 )5. 


ered in a scalar relativistic approximation. 

The enhancement of the SOC strength was per¬ 
formed for the ground state configuration of this 
system with lattice parameters a = 4.80 A and 
optimized c. The evolution of the gap at both k- 
points M and K with the ratio A so /Ao (Ao being the 
physical value of the SOC strength) is presented in 
Fig. [5] (b). The gap becomes smaller with increas¬ 
ing SOC strength and it closes and reopens again 
at A S o/Ao « 2.1. The topological invariants of the 
band structure Vb remain the same at both sides of 
the zero-gap point. As a consequence of this, the 
insulating phase remains trivial and no TI state was 
found. However, a discontinuous change in the gap 
occurs at that A so /Ao, as we can see in Fig. [5] (b). 
Further studies are required to understand this in¬ 
triguing behavior. 


D. Evolution of the gap with U at constant 
lattice parameters 

The 5 d electrons can exhibit different degrees of 
correlation effects. As explained above, the typical 
values of the on-site Coulomb repulsion can be on 


the order of 1 — 3 eV. Still, even if the actual U 
value is difficult to determine accurately, it is inter¬ 
esting to study the evolution of the band gap and 
the topological properties of the system as a func¬ 
tion of the Coulomb on-site repulsion. To take this 
effect into account, LDA+U calculations were per¬ 
formed for both fillings analysed above retaining a 
non-magnetic solution. 


1. LDA + U calculations for the {Au 2 0 i)/{Ah 03 )s 
system 

For the 5 d 8 configuration, the evolution of the gap 
with the Coulomb interaction was studied for three 
fixed in-plane lattice parameters: Qaz 2 0 3 = 4.89 A, 
a = 4.76 A and a = 4.95 A. First of all, even though 
now the calculations include an uncorrelated part 
treated with LDA, the results are consistent with the 
previously presented GGA calculations (Fig. [6] (a)), 
the topological phase of those calculations matches 
the ones obtained now at U = 0 using the LDA. The 
two first lattice parameters (see Fig. 11 were found 
to be topological insulators within GGA, and the 
system with a = 4.95 was characterised as a trivial 
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FIG. 6 . (Color online) (a) Evolution of the gap with the on-site Coulomb repulsion U for the system (AU 2 O 3 ) with 5d 8 
filling for the lattice parameters indicated. A transition between a topological insulator an a trivial one is produced 
at U = 2 — 3 eV for the studied configurations, (b) Band structure calculation for the Au 203 -based multilayer (5d 8 
filling) with a = 4.76 A for U = 2.31 and U = 2.72 eV. The increase in U induces a change in the topological 
properties of the system from a low-U topological phase to a higher-U trivial insulating solution. All of the LDA+U 
calculations performed show a non-magnetic solution. 


insulator, and so are now at zero U (using LDA). 

The evolution of the gap with the on-site Coulomb 
repulsion U is shown in Fig. [ 6 ] (a). The configura¬ 
tions with lattice parameters aAi 2 o 3 = 4.89 A and 
a = 4.76 A are topological insulators at U=0 and 
this topological nontriviality holds until their gaps 
close and reopen in the interval of 2 — 3 eV to a 
high-U value phase that is a trivial insulator. The 
configuration of a = 4.95 A is slightly different. At 
U = 0 it was found to be a trivial insulator, but in¬ 
creasing U the system evolves first into a topological 
insulator (very near U = 0) behaving in the region 
of 0.5 — 3 eV very similarly to the other two lat¬ 
tice parameters explored. The zero gap point, where 
the phase transition to a high-U trivial insulator oc¬ 
curs, is shifted to lower U values (from U = 2.45 to 
U = 2.18 eV), as tensile strain is applied within the 
range discussed here. 

Band structure calculations were performed for 
a = 4.76 A at both sides of the closing gap and 
these are shown in Fig. [ 6 ] (b) (the other two cases 
are analogous but presenting their phase transitions 
at different U values). By calculating the Vb of each 
band, we have identified that the system is a topo¬ 
logical insulator before the gap becomes zero and 


as U increases, the gap reopens and the system be¬ 
comes a trivial insulator at larger U. Introducing 
static correlations at the LDA+U level does not 
change drastically the electronic structure. Physi¬ 
cally meaningful Coulomb interaction values (1—3 
eV) host non-trivial topological phases and larger 
correlations tend to bring about the disappearance 
of the topological properties, as a general rule, quite 
strain-independently. 


2. LDA + U calculations for the (OS2O3)/(AfeChA 
system 

For the 5 d 5 filling, LDA+U calculations were per¬ 
formed for the lattice parameter of the AI 2 O 3 sub¬ 
strate CLAI 2 O 3 = 4.76 A and the ground state a = 4.80 
A. Like in the previous GGA calculations and SOC 
enhacement, no TI state was found. For both cases, 
the evolution of the gap with U (presented in Fig. |T| 
(a)) is quite similar. The large trivial gap closes and 
reopens again as a trivial insulating state in the re¬ 
gion of 4.7 — 5.3 eV. The way the gap does so affects 
both the M and K points in the band structure, 
leading to an intermediate metallic phase. 
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FIG. 7. (Color online) (a) Evolution of the gap versus the on-site Coulomb repulsion in the (OS 2 O 3 )/(AhC^s system. 
In contrast with previous studies of 5 d 5 filled honeycomb lattices, no configuration with a non-trivial topology was 
found with this corundum-based structure. The gap closing and reopening occurs by crossing through an intermediate 
metallic phase, (b) For the 5d 5 , the evolution of U produces a reordering of the topological invariants of the band 
structure filling at a = 4.80 A. All of the LDA+U calculations performed show a non-magnetic solution. 


In the ground state configuration, a reordering of 
the topological invariants of the f 2g manifold was 
found not for bands just around the Fermi level but 
for the lower lying t 2 g bands. We can see in Fig. Q] 
(b) the evolution of the ground state on both sides 
of the transition. A trivial phase is always obtained. 
However, the topological invariants of each band dif¬ 
fer, indicating there is a phase change induced by 
U, but both phases being trivial. For all physically 
meaningful U values, the system is in the same triv¬ 
ial insulating phase, the transition occurring at large 
unphysical correlation strength values. 


V. SUMMARY 

All along this work we have shown the conditions 
that allow the existence of non-trivial insulating 
phases for a metal-based honeycomb lattice created 
out of the corundum multilayered (M 2 03 )/(A 1 2 03)5 
system. This is a realization of an oxide-based hon¬ 
eycomb lattice of metal cations that does not require 
the use of unstable polar surfaces for its growth. 

We have studied different 5 d fillings that could 
present non-trivial phases by analyzing their band 
structure around the Fermi level, the existence of 


Dirac cones and the effects of adding spin-orbit cou¬ 
pling. A Dirac cone near the Fermi level was ob¬ 
tained for the Au-based (at 5 d 8 filling) and also for 
the Ir-based systems (shifting the chemical potential 
to one electron less at 5 d 5 filling), which open a gap 
with SOC. 

GGA calculations were performed in order to 
study the stability of the emerging non-trivial phases 
with in-plane strain. We observed that the 5 d 8 fill¬ 
ing is a TI with the lattice parameters of Al 2 03 , 
but with increasing a (tensile strain), the gap closes 
and reopens in a trivial way. The same procedure 
was employed in the 5d 5 filling but no TI phase was 
found. For this case (Os-based compound) we tested 
an enhancement of the SOC strength, and yet again 
all the systems calculated were found to be trivial 
insulators. 

On-site Coulomb repulsion effects for the 5d TM 
cations in the honeycomb were studied by means 
of the LDA+U method, for testing the appear¬ 
ance of magnetism and studying the phase diagram 
in terms of the topological properties of the sys¬ 
tem. Considering only non-magnetic solutions, the 
(Au 2 03 )/(A 1 2 03)5 compound presents a phase tran¬ 
sition from a low-U TI to a higher-U trivial insula¬ 
tor for values of U between 2 — 3 eV. However, the 
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( 0 s 2 C> 3 )/(Al 203)5 compound was found always in a 
trivial state, independent of the U value utilized. 

We have also studied the robustness of these ma¬ 
terials with respect to magnetic ordering. We found 
that AF solutions are eventually stable for both com¬ 
pounds but they appear at different values of U, our 
DFT calculations predict that, due to the small band 
widths of the e g states, an AF order would be present 
for every realistic U value 1 — 3 eV in the 5 d 8 con¬ 
figuration, whereas in the 5 d 5 system much larger U 
values (higher than 5 eV) are necessary to stabilise 
it. The existence of a magnetic phase based on Au 3+ 
cations would be rather unique, caused by the band 
width reduction due to spatial confinement in those 
nanostructures. 

The results for the time-reversal symmetric phase 
at 5 d 8 filling in terms of topological properties are 
similar to those found in buckled honeycomb lat¬ 
tices created out of perovskite bilayers, where non¬ 
trivial phases were also predicted. On the other side, 


the 5 d 5 filling behaves differently, and no topologi¬ 
cal phases appear. Also the band structures are dif¬ 
ferent, with the Dirac cones appearing at different 
places in the Brillouin zone. Thus, we can conclude 
that although in both cases Dirac cones appear at 
the Fermi level, the different hopping paths in the 
corundum structure compared to the perovskite bi¬ 
layers yield a different tight binding model with dis¬ 
tinct topological phases. 
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